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ABSTRACT 


A generalization  of  Freund's  bivariate  exponential  model 
is  discussed.  Probabilistic  properties  of  this  model  with 
minimal  distribution  assumptions  are  derived  including  the 
joint  survival  function  and  Laplace-Stielt j es  transform. 
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1.  Introduction.  A study  is  made  of  a two  component 
failure  system.  Let  the  components  be  denoted  by  A and  B 


and  their  lifetimes  by  S and  T respectively.  The  general  model 
to  be  discussed  in  this  paper  is  a generalization  of  that  of 
Freund  |2]  and  arose  from  consideration  of  bivariate  exponential 
distributions.  In  addition  to  Freund's  there  have  been  many 
such  distributions  suggested  in  the  literature,  including  those 
suggested  by  Gumbel  [3],  Downton  [1]  and  Hawkes  [4].  Marshall 
and  Olkin  [5]  have  proposed  the  most  widely  referenced  model. 

In  their  model,  the  Pr  (S>s+A | S>s ,T>s } = Pr  {S>s+A | S>s } . This 
implies  that,  conditioned  on  the  fact  that  component  A is 
functioning  at  time  s,  the  distribution  of  its  residual  life- 
time is  independent  of  whether  component  B has  failed  or  not. 

Freund  [2]  derives  his  distribution  from  the  assumption 
that  at  the  failure  of  one  component  the  distributio " of  the 
residual  lifetime  of  the  other  component  is  changed.  Specifically, 
he  assumes  that  the  lifetime  of  A given  that  no  failure  of  B 
occurs  is  distributed  exponentially  with  parameter  a,  denoted 
exp  (a),  and  that  the  lifetime  of  B given  no  failure  of  A is 
distributed  as  exp  (B).  If  B fails  before  A,  then  the  residual 
lifetime  of  A is  distributed  as  exp  (a').  If  A fails  before 
B,  then  the  residual  lifetime  of  B is  distributed  as  exp  (8'). 

This  can  be  summarized  as  follows.  Let  X,Y,U,V  be  independent 
random  variables  such  that  X ~ exp  (a),  Y ~ exp  (B) , U - exp  (a') 
and  V ~ exp  (8').  Then 


2 


( x 

if 

X < Y 

[x+v 

if 

X 

< Y 

s = \ 

~ , T - 

1 

1 

[Y+U 

if 

X > Y 

[ Y 

if 

X 

> Y 

The  model  introduced  in  section  2 is  a generalization 
of  Freund's  model  with  the  distribution  assumptions  dropped. 
Section  3 includes  derivations  of  expressions  for  the  joint 
survival  function  and  marginal  distributions  for  this  model. 
Letting  Y ~ exp  (8),  an  explicit  expression  for  the  joint 
Laplace-Stielt jes  transform  of  the  distribtuion  is  found  in 
section  4.  Using  this  express  ion,  moments  of  the  distribution 
are  calculated.  In  section  5,  two  definitions  of  system  life 
are  defined  and  the  Laplace-Stielt j es  transform  of  the  system 
life  distribution  is  found  in  each  case. 

2.  Model  Definition.  Label  the  two  components  of  a 
system  A and  B.  The  lifetimes  of  the  two  components  are 
dependent,  in  that  the  failure  of  one  component  affects 
the  residual  lifetime  of  the  other.  Formally,  let  A and  B 
have  lifetimes  S and  T respectively.  Given  random  variables 
X,Y,U,V  write 
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to  represent  the  two  lifetimes.  The  following  assumptions  are 
made  on  X,Y,U,V: 
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i)  X,Y,U,V  are  mutually  independent, 

ii)  Pr  { X >0}  = 1,  Pr  {Y  >0}  = 1, 

iii)  Pr {U  <0}  = 0,  Pr(V  <0}  = 0, 

iv)  X,Y  are  absolutely  continuous. 


(2.2) 


Let  Pr(U  = 0}  be  denoted  by  Py  and  Pr{V  = 0}  by  pv . The 
following  is  an  example  of  the  applicability  of  the  model. 

Let  A and  B be  two  generators  that  supply  electricity  to 
a hospital,  each  supplying  a portion  of  the  building.  When 
one  generator  fails,  the  remaining  gererator  must  supply  power 
to  the  entire  hospital.  This  puts  added  strain  on  the  generator 
and  thus  alters  its  residual  life. 

3.  The  Joint  Survival  and  Marginal  Distributions.  Let 
F(s,t)  = Pr(S>s,T>t}.  We  begin  by  deriving  an  expression  for 
F(s,t)  when  s < t.  First  condition  on  the  values  of  X and  Y. 

.00  oo 

Then  F(s,t)  = Pr(S>s,  T>t  | X=x,  Y=y}  dFY(x)  dFY(y) . Now 

J0  Jo  X 

partition  the  region  [0,®)  * [0,®)  into  12  subregions  based  on 
the  relative  sizes  of  x,y,s,t.  The  only  non-zero  contributions 
are  Pr{V>t-x}  on  the  region  [s,t]  x [x,®)  and  1 on  the  region 
[t,®)  x [t,®) . Therefore  when  s < t, 

>X  f°°  f°°  f°° 

F(s,t)  = Fy(t-x)  dFy(y)  dFx(x)  + 1 dFx(x)  dFY(y) 


Fv(t)F  (t)  + F.r(t-x)  Fy(x)  dFY(x) 


L 


The  case  s > t follows  similarly  and  the  case  s = t is  trivial 
We  have  therefore  proved 


Theorem  1 


We  can  obtain  the  marginal  distributions  using  Fg(s)  = F(s,0) 


Corollary  2 


From  the  model  definition  we  have  the  following  representation  for 


the  variable  S 


S = min(X.Y)  + U • I 


where  ^ { x > Y } *s  t*ie  inc*icator  function  for  the  set  {X>Y}.  A 
similar  representation  exists  for  T.  With  this  representation  we 
see  that  in  general 


E (S)  = E (min  (X  ,Y) ) + E(U)  Pr{X>Y} 
E (T)  = E (min  (X , Y) ) + E(V)  Pr{Y>X} 


4.  The  Joint  Laplace-Stielt jes  Transform.  The  joint  Laplace 


Stieltjes  (L-S)  transform  of  (S,T)  is  defined  to  be 

r r -as-bt 


F(ds,dt),  where  F(ds,dt)  is  the  measure  determined 
by  the  survival  function  F(s,t).  We  begin  with  a representation 


Lemma 


Proof  for  the  case  s<  t:  We  must  investigate  how  A can  fail  at 


time  s and  B fail  at  time  t.  If  set,  then  S < T,  and  thus  X<Y 


so  that  S = X = s and  T = X + V 


and  V=t  = s.  Since  the  variables  are  independent  the  result 


follows 


Using  this  lemma,  the  joint  L-S  transform  of  (S,T)  can 


be  written  as 


This  expression  can  be  evaluated  piece  by  piece.  We  have 


- as  - d 
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Here  f*(b)  is  the  L-  S transform  of  Fy  evaluated  at  b. 


By  a similar  calculation 


-as-bt  tt 


0 Jt  + 


e Fx(t)  dFy(s-t)  dFy(t) 


= (ffl  (a)  - Py) 
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Fv(s)  dFv(s)  . 
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Combining  these  into  (4.1)  we  get 
Theorem  4. 


e (a+b)s  Fx(s)  dF y( 


f*(a,b)  = fy(b)  ( e b-*S  Fy(s)  dFx(s)+  fy(a) 

J0 

These  integrals  cannot  be  evaluated  in  general,  but  can  be  for  certain 
important  special  cases.  In  particular  we  have 

Corollary  5. 

If  Y ~ exp  (g>  then 


f*(a)6 


+ a+b  + B 


f*(a,b)  = ffc(b)  f$(a+b-*6  ) 


[ 1 - f*  (a+b  + B) ] • 

A 
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Corollary  6. 

If  X ' exp  (a)  and  Y ~ exp  (6),  then 

f*(a,b)  - [«f{(b)*6fj(a)]. 

Using  the  expression  in  Corollary  6,  moments  of  (S,T)  can  be  found 
by  evaluating  the  appropriate  partial  derivatives. 

Theorem  7. 

If  X ~ exp  (a)  and  Y-exp  (6),  then 

i)  E (S)  = ^ [l  + |SE(U)]f  E(T)  = ^ [1+aE  (V) ] , 

ii)  Var(S)  - 1 , U‘82  Var(U).ceE(U2))  , 

(a+8) 

Var(T)  = — *•  [ 1 + a2  Var (V) +aBE (V2) ] , 

(a+B) 2 

iii)  Cov(S,T)  = - — j [l-aBE(U)  E (V)  ] . 

(a+B) L 

It  can  be  shown  that  all  values  between  +1  and  -1  are  possible  for 
the  correlation  of  S and  T even  in  the  special  case  considered  in 
Theorem  7. 

5.  Time  to  System  Failure.  Failure  in  a series  system  and  a parallel 
system  will  be  discussed.  Label  the  system  life  L.^  for  the  series 
system  and  L2  for  the  parallel  system.  We  have  = min  (S,T)  and 
L2  = max  (S,T) . 

Since  min  (S,T)  = min  (X,Y) , F,  (t)  = F ft)  F ft)  and  the  model 

l j x y 

does  not  really  enter  in.  If  X and  Y are  exponentially  distributed 
then  so  is  L,. 
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To  work  with  I^,  the  random  variable  D = will  be  considered. 

If  we  let  p = Pr  (X  > Y}  and  q = 1 - p,  then  D ~ pV  + qU.  Since  V and 
U are  independent  f*  (a)  = pq  f*  (ap)  f*  (aq)  • It:  follows  from 
L2  = + D and  the  independence  of  and  D that 

f*L  (a)  = Pqf * Cap)  f*(aq)  f *L  (*)  • 

6.  Discussion.  The  question  of  practicality  arises  here  as  it 
should  with  any  model  formulation.  Is  the  class  of  models  sufficiently 
rich  yet  simple  enough  to  obtain  meaningful  results?  The  proposed 
model  is  based  upon  the  assumption  that  the  residual  life  of  one 
component  is  dependent  upon  whether  or  not  the  other  component  has 
failed.  This  seems  to  be  a realistic  assumption  in  many  applications. 

The  first  meaningful  results  one  would  like  are  the  moments  of 
(S,T).  For  given  distributions  the  Laplace- Stieltj es  transform 
(Theorem  4)  may  be  used  to  obtain  the  moments.  This  expression  may 
be  difficult  to  evaluate  though.  However,  the  much  simpler  expression 
in  Corollaries  5 and  6 still  come  from  rich  classes  of  models. 

Next,  given  a sample  can  one  obtain  reasonable  estimators  for  the 
parameters?  In  Freund's  special  case,  simple  expressions  for  the 
maximum  likelihood  estimates  can  be  obtained. 

Finally,  does  this  model  enable  one  to  make  decisions  about  how 


to  control  the  system.  When  the  components  fail,  can  we  arrive  at 
optimal  replacement  policies?  Here  again,  this  is  possible  when 
using  this  model. 
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